Influence of thermal coupling on spin avalanches in Mni2-acetate 
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The effect of thermal couphng on spin avalanches in Mni2-acetate has been probed using a single 
crystal assembly. Time-resolved, synchronized measurements of magnetization and temperature are 
reported. Unusually low avalanche trigger fields occur when thermal coupling to the bath is weak. A 
temperature rise observed at zero magnetic field is attributed to a change in magnetostatic energy. 
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Over the past decade there has been much interest in 
single molecule magnets due to their macroscopic quan- 
tum behaviour. Mni20i2(CH3COO)i6(H20)4, abbrevi- 
ated to Mni2-acetate or Mni2, is one of the most stud- 
ied single molecule magnets to date.^ Strong uniaxial 
anisotropy and a large spin [S = 10) enable macroscopic 
quantum tunneling (MQT), in which the spin can reverse 
its direction by tunneling through the anisotropy bar- 
rier. In bulk crystals of Mni2 such tunneling events can 
be observed as steps in the magnetization curve, which 
occur at regular intervals in the external magnetic field 
Bn = uBq, where n is an integer and Bq = 0.46 
Recently, however, much attention has focussed on spin 
avalanches, in which the entire ensemble of spins reverses 
very rapidly.^ Spin avalanches are generally accepted to 
be thermal phenomena, driven by the cooperative pro- 
cess of spin-phonon relaxationi^i They occur only when 
the sample is sufficiently large that phonons emitted 
by the reversing spins cannot escape to the bath. The 
most recent research has shown that avalanches propa- 
gate through the crystal at constant velocity — a classical 
phenomenon dubbed "magnetic deflagration".^ 

In this Brief Report wc study the influence of ther- 
mal coupling to a heat bath on spin avalanches in a 
single crystal assembly of Mni2. We report measure- 
ments of the magnetic hysteresis curve by superconduct- 
ing quantum interference device (SQUID) magnetome- 
try and time-resolved, synchronized measurements of the 
sample temperature and magnetization, using a pickup 
coil. The effect of thermal coupling on the trigger field 
is discussed. We also study a temperature rise at zero 
magnetic field that is not associated with avalanches. 
Similar temperature rises have been observed by other 
researchers but have never been satisfactorily ex- 

plained. 

The experiments described in this Brief Report were 
performed on an assembly of single crystals of Mni2- 
acetate grown according to the procedure described by 
LiSf^ Each crystal had a length ~ 2 mm and a cross 
sectional area ~ 0.5 x 0.5 mm^. The long dimension is 
known to correspond to the c-axis. Hence, we were able 
to glue 35 crystals together using G. E. varnish to create 



an assembly of mass to = 21 mg with the c-axes aligned 
to within ~ 5°. 

We began by measuring two sets of magnetic hystere- 
sis curves using a commercial SQUID magnetometer. In 
the first set, the sample was strongly coupled to the bath 
via a constant flow of helium exchange gas. In the sec- 
ond set, the thermal coupling was weakened by enclosing 
the sample in an evacuated glass tube. For both sets 
the sample was mounted with the c-axis parallel to the 
external magnetic field to within 5°. 

For the first set of measurements we initially cooled 
the sample to 1.8 K with the external magnetic field set 
to zero. By ramping the field slowly to 2 T we found 
the saturation moment^^ to be 0.95 mJT"-'^. We then 
measured a series of hysteresis curves at different tem- 
peratures, each time sweeping from +2 T to —2 T and 
back at a rate of 3 mTs~^. Avalanches were consistently 
observed for T < 2.4 K, but the field values at which they 
were triggered were not completely reproducible, giving 
rise to asymmetric hysteresis curves [Fig. HJa)]. In ad- 
dition, the trigger fields did not always correspond to 
the resonance fields B^. This is in contrast to previous 
measurements made on similar-sized assemblies of single 
crystals at similar temperatures, ^'^'^^ but in agreement 
with measurements made on individual single crystals at 
lower temperatures.— Above 2.4 K, the avalanches disap- 
peared and MQT steps were recovered. Above 3.6 K, the 
magnetization curve became non-hysteretic. 

For the second set of measurements the thermal cou- 
pling was weakened by enclosing the sample in an evac- 
uated glass tube that functioned as a shield against the 
flow of helium exchange gas. It was necessary to reduce 
the sample volume by ~ 20% in order to flt it inside the 
tube. This reduced the mass from 21 to 17 mg and the 
saturation moment from 0.95 to 0.77 mJT^^. The rest 
of the experiments described in this Brief Report were 
conducted on the diminished sample. 

Again, we measured hysteresis curves at a series of 
temperatures [Fig.[TJb)]. As before, spin avalanches were 
observed for T < 2.4 K. However, this time the trig- 
ger fields were lower and more reproducible, leading to 
symmetric hysteresis curves. In addition, the magnetic 
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FIG. 1: (Color online). Magnetic hysteresis curves measured 
by SQUID magnetometry. The magnetic moment is normal- 
ized to its saturation value /is at B = 2 T. (a) Strong thermal 
coupling between sample and bath. Avalanches are observed 
between 1.8 and 2.4 K. Resonant MQT steps are observed 
between 2.7 and 3.3 K. (b) Weak thermal coupling between 
sample and bath. Avalanches are observed between 1.8 and 
2.4 K. Above 2.4 K the magnetization curve is smoothly hys- 
teretic with no discontinuities. 



moment approached saturation more slowly leading to 
rounding of the avalanche steps as they approached sat- 
uration. Above 2.4 K no MQT was observed, in contrast 
to the multiple steps observed with the strongly coupled 
sample. Instead, the hysteresis curve became smoothly 
hysteretic. The loss of MQT in the weakly coupled sam- 
ple suggests a strong rise in sample temperature due to 
spin-phonon relaxation. When thermal coupling to the 
bath is strong, heat is rapidly dissipated, but when the 
coupling is weak the heat cannot escape so quickly to the 
bath. It appears that the glass tube provides sufficient 
thermal isolation that sample heating causes the spin re- 
versal mechanism to be dominated by thermal activation, 
rather than by MQT. 

There are several reports in the literature of sig- 
nificant rises in the sample temperature during spin 
avalanches^ii^iii^ and MQT)2iiiSiii with reported values 
ranging between 1 and 5 K. The spread is probably due 



to variations in sample size, field-sweep-rate, location of 
the thermometer, and thermal coupling between the sam- 
ple and the bath. Here, we report time-resolved mea- 
surements of the temperature rise during avalanches in 
a weakly coupled sample. Our measurements were per- 
formed in a helium cryostat with a superconducting mag- 
net. We placed the sample at the bottom of an open- 
ended glass tube that functioned as a weak thermal link 
to the bath. A multiturn coil was wound around the 
base of the tube to enable us to detect changes in the 
sample magnetization. The tube was then placed inside 
an evacuated copper sample holder with a vertical slit to 
minimize eddy currents, which was attached to the 1 K 
plate. We estimate that the c-axis of the crystal assem- 
bly was aligned parallel to the magnetic field to within 
~ 10°. A Cernox resistance thermometer was mounted 
directly on top of the sample using high thermal conduc- 
tivity grease ( Apiczon N) . A low bias current of 5 nA was 
used to minimize self heating. 

The sample was initially cooled to 1.6 K with the ex- 
ternal magnetic field set to zero. We then swept the 
field at a rate of 24 mTs^^ and used a VXI data aqui- 
sition board to record time-resolved signals from either 
the coil or the thermometer. Separate measurements 
were subsequently synchronized by aligning their trigger 
points. The temperature and magnetization measure- 
ments were both triggered from the thermometer signal, 
which was found to be highly reproducible on successive 
field sweeps. The rate of change of the sample magneti- 
zation dM/dt was calibrated by measuring the constant 
voltage induced across the coil by the external field ramp 
dB/dt. This was done at high fields (between 2 T and 
3.5 T) where the sample magnetization was saturated. 

Figure [2I^a) shows the changes in sample temperature 
and magnetization observed during an avalanche. The 
oscillation in the temperature signal is due to electri- 
cal pickup of the 50 Hz line voltage A distribution of 
avalanche trigger fields was observed between 0.37 and 
0.65 T, with a mean of 0.56 ± 0.09 T, which is close to 
Bi. The temperature rise was found to be quite repro- 
ducible, despite the spread in trigger fields, with a mean 
maximum temperature of 4.3 ± 0.1 K. 

In order to compare the time-resolved avalanche mea- 
surements with our earlier SQUID measurements we in- 
tegrated the calibrated signal dM/dt and multiplied it 
by the sample volumei^ to obtain the change in mag- 
netic moment Afj,. This yielded A/j, = 0.03 mJT~^ for 
all avalanches, regardless of the trigger field, which is a 
great deal smaller than the change in moment A/i « 1.5 
mJT~^ observed in the SQUID measurements on the 
weakly coupled sample. We believe the discrepancy to 
be due to the difference in timescale of the two measure- 
ments. As can be seen in Fig. [H^a), the spin reversal 
rate declines rapidly as the temperature rises, resulting 
in the termination of the avalanche after a few millisec- 
onds. However, the SQUID measurements on the weakly 
coupled sample show that spin reversal continues over 
a much longer period of ^ 100 s at a rate dM/dt < 2 
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FIG. 2: (Color online). Time-resolved, synchronized mea- 
surements of magnetization and temperature on a sample 
weakly coupled to the bath, (a) Main plots: spin avalanche 
at S = 0.65 T. The grey line is a guide to the eye. Upper 
inset: separate avalanche recorded over a longer time interval. 
Prior to this avalanche, a small rise in temperature occurs at 
B ~ Bi. Lower inset: distribution of trigger fields, (b) Step 
at B = 0. This plot comprises four separate measurements 
(two of temperature and two of magnetization) which have 
been aligned by their trigger points. Inset: small change in 
magnetization at B Ri B4,. 



kJ m~'^ (close to the noise floor of the coil). This 
suggests that the majority of the spin reversal detected 
by the SQUID was due to thermal activation. 

Besides avalanches, we also observed slower changes in 
the temperature and magnetization at i? = [Fig. [2Ub)]. 
The change in magnetization at _B = is known to arise 
from thermally activated relaxation of a minority molec- 
ular species constituting about 5% of the sample. The 
relaxation rate of this species is much faster than that 
of the majority species because the anisotropy barrier is 
lower. However, the rate of change of the total magnetic 
moment is still slow compared to an avalanche, because 
thermal activation is a non-cooperative process. Over 
a period of 20 s, the change in magnetic moment was 
found to be A/i = 0.06 mJT~^, in good agreement with 
the SQUID measurements over a similar timescale. The 



temperature rise was found to be very reproducible, with 
a mean maximum temperature of 2.2 ± 0.05 K. 

In addition to the above observations, we also ob- 
served occasional small signals before or after the main 
avalanche. Small temperature rises of ^ 50 niK were ob- 
served a.t B Bi [insets in Figs.[2lja) and[2ljb)] and small 
changes in magnetization were observed close to B « B3 
and B4. Because these signals were observed only occa- 
sionally, we were unable to synchronize temperature and 
magnetization measurements. Therefore, we cannot tell 
whether they were due to MQT or miniature avalanches 
involving parts of the sample in poor thermal contact 
with the rest. However, the duration of the magnetiza- 
tion signals 100 ms) suggests that they are due to 
MQT rather than miniature avalanches, which in turn 
suggests that the main avalanche does not always result 
in full spin reversal. 

We now report the results of a second time-resolved 
experiment, this time with the sample enclosed in a glass 
tube filled with helium exchange gas to improve the ther- 
mal coupling. The experimental conditions were identi- 
cal to those of the previous experiment, except that we 
were unable to measure the sample temperature, due to 
the difficulty of sealing the glass tube around the sensor 
leads. Avalanches were observed at a variety of trigger 
fields between 0.91 and 1.42 T [Fig. [Sla)] . These corre- 
spond approximately to B2 and B3. The trigger field was 
consistently higher than when the glass tube was evacu- 
ated [Fig. [3^b)], in agreement with our SQUID measure- 
ments. However, the change in magnetic moment A/i 
was similar in both time-resolved experiments, suggest- 
ing that, in both cases, avalanches terminate long before 
full spin reversal is achieved. This is supported by the 
observation of two closely spaced bursts of spin reversal 
at the lowest trigger fields in the evacuated sample [inset 
in Fig.Etb)]. 

We determined the avalanche duration by fitting a 
Gaussian profile and taking the full width at half 
maximum. The duration was found to decrease expo- 
nentially as the trigger field increased, independent of 
thermal coupling to the bath [inset in Fig.[3l^a)]. This is 
consistent with the recently proposed model of magnetic 
deflagrationj^ii^ A particularly interesting observation is 
that, in the evacuated sample, avalanches triggered at 
identical fields had identical spin-reversal patterns. Since 
the lack of smoothness is probably due to sample inho- 
mogeneity resulting from the layers of varnish that bind 
the assembly, this suggests that avalanches propagate in 
a very reproducible manner. 

In the final part of our report, we investigate the cause 
of the temperature rise at _B = 0. The amount of heat 
generated is given by AQ = mcAT, where m is the mass 
of the sample, c is the specific heat capacity and AT 
is the observed temperature rise. Fominaya et alJ^ have 
measured the specific heat capacity of Mni2 as a function 
of temperature at zero magnetic field. Using a polyno- 
mial fit to extrapolate their data to temperatures below 
3 K we obtain c = 0.8 Jkg-^ K'^ at T = 1.6 K (the ini- 
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FIG. 3: (Color online), (a) Avalanches in sample with helium 
exchange gas (offset vertically for clarity). Inset: avalanche 
duration vs trigger field. Open circles: sample with exchange 
gas; closed circles: evacuated sample; solid line: exponential 
fit. (b) Avalanches in evacuated sample. Inset: two bursts of 
spin reversal are seen at the lowest trigger fields. 



tial temperature) and c = 1.3 Jkg~^K~^ at T = 2.2 K 
(the peak temperature). Taking the average of these, we 
obtain (c) « 1.1 Jkg~^K~'^. The observed temperature 
rise was AT = 0.6 K and the sample mass was m — 17 



mg, leading us to estimate AQ « 11 fiJ. 

At B = 0, no work is done by the external magnetic 
field because J B d/i = 0. However, work can be done by 
the internal demagnetizing field _Bd = —noNfi/V, where 
iV is the demagnetizing factor, /xq is the permeability 
of free space and V is the volume of the sample. The 
magnetostatic energy Em = J B^dfi is given by 



Em — 



2V 



(1) 



Taking the sample to be approximately spherical we as- 
sume a demagnetizing factor N = 1/3. The total change 
in the magnetic moment between full saturation and 
the beginning of the avalanche is Afj, 0.2 mJT~^. 
The amount of magnetostatic energy released is there- 
fore AEjn w 10 fj,3, in good agreement with AQ. 

In conclusion, we have compared avalanche measure- 
ments in samples with weak and strong thermal coupling 
to a heat bath and have shown that a decrease in thermal 
coupling leads to a reduction in trigger field. This sup- 
ports the view that heating is an important factor in the 
triggering of avalanches. Time-resolved measurements of 
the sample magnetization have shown an exponential de- 
pendence of the avalanche duration on the trigger field, 
consistent with the model of magnetic deflagration. A 
comparison of time-resolved measurements with SQUID 
magnetometry measurements has shown that avalanches 
terminate long before saturation is reached. The ele- 
vated sample temperature following an avalanche allows 
spin reversal to continue via thermal activation. 

In addition to our studies of spin avalanches, we have 
also observed significant heating at _B = 0. This has 
been observed previously by other researchers,i^'iS,!il but 
no satisfactory explanation has been given. We attribute 
the temperature rise to a release of magnetostatic energy. 
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Throughout this Brief Report, we use the symbol fi for the 



magnetic moment of the sample and Af for the magneti- 
zation (M = IJ./V , where V is the sample volume). 
The amplitude of the 50 Hz oscillation gets larger as the 
temperature increases because the sensor resistance de- 
creases by an order of magnitude between 1.6 and 4.3 K. 
This causes the ratio of pickup voltage to sensor voltage to 
decrease. 

The sample volume = 5.5 mm^ was computed from the 
saturation moment fj.3 = 0.77 mj T~^ measured by SQUID 
magnetometry, using the known crystal structure and unit 
cell volumsi^ and the known spin S = 10 per molecule. 
In the case of the double peaks observed at the lowest 
trigger fields, only the first peak was fitted. 



